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We establish a comparison between Rakib-Sivashinsky and Michelson-Sivashinsky quasi- 
linear parabolic differential equations governing the weak thermal limit of upward flame front 



^ . propagating in a channel. For the former equation, we give a complete description of all steady 

H ; 

5t , solutions and present their local and global stability analysis. For the latter, multi-coalescent 

unstable steady solutions are introduced and shown to be exponentially more numerous than 
the previous known coalescente solutions. This fact is argued to be responsible for the dis- 
agreement of the observed dynamics in numerical experiments with the exact (linear) stabil- 
ity analysis and also gives the ingredients to describe the quasi-stable behavior of parabolic 
steadily propagating flame with centered tip. 
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Rakib-Sivashinsky quasilinear parabolic differential equation governing the weak thermal limit 
of an upward flame interface propagating in a channel is given by 

(f>t = e(l)xx - ^(pl + 4>- (i>, (1) 
where y = (j){t,x), < x < tc, defines an instantaneous flame profile in dimensionless variables, 

- 1 r 

0= — / (j)(t, x) dx denotes the space average, e > is a physical parameter (Markstein length) 
Jo 

and Neumann (adiabatic) boundary condition at channel walls is imposed: (f)x{t,0) = vr) = 0. 
According to Darrieus-Landau's hydrodynamic flame theory |]R^] there are two competing sources 
of instabilities in this model given by the last two terms of (|l|) . A plane flame front separating the 
cold and hot gas is subject to Rayleigh-Taylor instabihty, due to the thermal expansion, and to the 
buoyancy effect caused by external acceleration. These two ingredients lead the front to become 
convex toward the cold gas. As in |[RS|| , the acceleration vector points toward positive y-values and 



the upward propagating direction is negative. For the purpose of comparison such convention will 
be maintained. 

Paraboloid profiles with the tip located around the center of the channel are observed, both 



experimentally and in computational simulations, to remain "stable" for long time (see e.g. ||MiS 
and references therein). There are, however, experiments whose paraboloid tip eventually slides to 
the channel wall. A mathematical description of the former configuration as a "quasi-equilibrium 



transient state" was provided by Berestycki, Kamin and Sivashinsky | [BKS| | in which work stationary 



solutions of (|I|), their respective stability properties and the nonlinear dynamics were described for 
sufficiently small e. Metastable dynamics has been studied in details in |^W|| . The proofs of the 
other theorems stated in pKS|| remain, to the best of our knowledge, unpublished. 



Sivashinsky's previous equation of weak thermal expansion[g, 

(f>t = - + I{4>) , (2) 

hasn't considered buoyancy effect and Darrieus-Landau instability has been taken into account 
replacing 0— in ([l|) by a linear singular integral operator given by a multiplication by \k\ in 
the Fourier representation 

I {(j)) {t,x) = - ^k / cos[k{x - y)] (t){t,y)dy . 
^ k=i -'-■^ 



In contradistinction, the only trace of term in equation (P comes from the removal of = 
mode 

4>-4>=-^ cos[k{x -y)] (t>{t,y)dy . 
^ k=i -^-^ 

Michelson-Sivashinsky equation (0) has been studied by many authors after Thual-Frish-Henon's 
application ||TFH| of pole decomposition. For a more recent survey of this method, see Vaymblat 



and Matalon |[VMl| , |VM2|| and references therein. 



Despite of the fact that equations (0) and (|^) differ in many respects, it is our purpose to expose 
the similarities and distinctions of their solutions. We shall see, by the introduction of so called bi- 
coalescent steady states, that equation may also admit quasi-equilibrium transient "parabolic" 
profile with centered tip. Our comparison relies on two recent analysis. Firstly, Vaymblat and 
Matalon [ |VM1|| have determined all steady solutions whose poles coalesce into a line parallel to the 



imaginary axis and solved explicitly the eigenvalue problem of (0) linearized about those states. 
They conclude, in a second paper IVM^] , the existence of only one linearly stable steady coalescent 
pole solution. 

For equation (|Tp with a different parametrization, we have given | |GM[ a complete description 
of all equilibrium solutions and provided their local and global stability in an appropriated Sobolev 



space. Since our results are major extensions of those stated in ||BKS|| and the mathematical 



presentation in ||GM|| may cause certain difficulties in translating to the present application, we 



shall here restate them with a brief explanations of their proofs. A detailed presentation including 
the metastable states analysis will appear elsewhere ||MG| . 



We shall first restrict ourselves to the Rakib-Sivashinsky equation. The details presented after 
each statement (detached by brackets) are essentials to establish the subsequent comparison but 
may be skipped in a first reading. Multi-coalescent steady states of Michelson-Sivashinsky equation 
will be considered next. A conclusion will be presented at the end. 

To discuss our results let us consider the equation 

Ut = -UUj: + U , (3) 

for the derivative u = (p^, with Dirichlet boundary conditions u{t,0) = u{t,'K) = 0. There is a 



one-to-one correspondence between 

eit,x) := / u{t,y)dy, (4) 







with u a solution of (||), and a solution of (|1|) given as follows. Notice that 9 satisfy (|T]) with 
replaced by e(j)xx(t,0). By definition, 9(t,0) = for alH > and their steady solutions do not 
propagate (they are equilibrium). We have 9(t,x) = 4>(t,x) — 0) , 9^ = (p^ and 9^^ = 'Pxx- For 
the opposite relation, we have 

(p+ [ ^{s)e'-'ds = 9 + e [ 9xx{s,0) e'-' ds , 
Jo Jo 

because both sides satisfy equation (|1|) without term, which yields 

<f){t, x) = 9{t, x)- {9 (r) - 6 9xx{t, 0)) ds . (5) 

This and @ will be used to discuss the steadily propagating solutions of (|l|). 

If A denotes the operator given by r.h.s. of (|^) linearized about the trivial solution uq = 0, let 

2 

be the Banach space equipped with the graph norm ||/||i/2 — / '^^ (equivalent 

Jo 

to the Sobolev space Hq{0,tt)). In ||GM|| , Theorems 3.2 and 5.10, we have proven: 

The initial value problem ^ with u{0, ■) = Uq E B^^"^ has a unique solution for all t > and the 
trajectories {u{t, ■)}^>o He in a compact set in B^^"^ . 

By @) and (|^), this provides the existence and uniqueness of the initial value problem ([l|) for 
alH > 0. Compactness property will be useful for the global stability analysis. 

The steadily propagating solutions of (|1]) can be obtained from the equilibrium solution of (|^). 
The quantity inside the parenthesis in (^ remains positive under the dynamics of (^ and the 
front propagates upward (toward negative y-values). As a consequence, the flame front propagates 
steadily with velocity 

V = - ^ +£?9x.(0) = -^^, (6) 
provided 'd[x) = v{y) dy, where f is a non trivial solution of 

ev^x - vv^ + V = , (7) 



with f (0) = f (vr) = 0. The next result is Theorem 4.1 of [pM 



For e > 1, vq = is the unique solution of (Q). For e < 1 such that ^ ^ ^ ^ £ < holds for some 
k & N, there exist 2k nontrivial solutions uf , j = 1, . . . , k, with j — I zeros in (0, vr) and such that 
u~{x) = —u^{tt — x) for j odd and u~{x) = u^(7r/j + x) mod vr for j even. Each pair uf bifurcate 
from the trivial solution vq = at Sj = with limu^ = 0. In the phase space M^, {u^' ,u^^ 
describes closed orbits around (0, 0) whose distance from the origin increases monotonically as e 
decreases. 

[Equation (|^ can be written as a dynamical system 

(8) 




where p = v and w = v. Since (^) remains unaltered by changing p — >■ —p and x — >■ —x, the orbits 
are symmetric with respect to the w-axis and Dirichlet boundary conditions p{0) = p{tt) = hold 
for any periodic orbit with period 2it such that p = aX x = 0. In ||GM|| we have shown that (i) 
is foliated by non-overlapping orbits 

= {{w{x),p{x)) : x G Mand (w(0),p(0)) = {wo,0)} 

labeled by the coordinate wq of the positive w-axis; (ii) as e varies, the orbit varies continuously 
from one to another; (ill) the orbits are closed if t^o < 1 and open if wq > 1; (iv) the origin (0, 0) 
is enclosed by any closed orbits. The trajectories 7„,(,, Wq > 0, are portrayed in Figure |l|. 

The period T = T{e, Wq) of a closed orbit varies continuously in e and Wq and is given by 

dp 



T = dx = 2 
'7^0 



w 



where the second integral is calculated in the semi-orbit. Hence, T = 2tx defines implicitly a 
continuous function wq = Wo{e), whose corresponding semi-orbit 7^^ = {{w{x),p{x)) : < x < vr} 
determines the solution v{x) = p{x) of (|^). The question to be addressed now is how many 2n- 
periodic orbits there are for each e. If Wq is small, 7^^ can be approximated by an ellipsoid: 



dv 

Tellipse = 4 / , ^ ^ ^ = 2n^ (9) 
Jo \/l - V 
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Figure 1: Trajectories of the dynamical system (^. 
uniformly in wq and we have 

lim T(5,wo) = 27rv^. (10) 

In ||GM|| , we have shown 

dT , , 

T^>0 (11) 

holds form any £ > 0. Consequently, T = T(e,wo) is a (strictly) monotonic increasing function 
of and e, since T(£:, wq) = a/£T(1,wo). As T increases when e and increase, equation ( p!OD 
implies that there is no nontrivial solutions (f 7^ 0) for £ > 1, since the period of each nontrivial 
orbit exceeds 27r. Now, let 

Gi(e, Wo) = r(e, wq) - 27r , 

be defined for < e < 1 and < < 1- Under the condition (|llD, there exist a unique function 
gx : [0, 1] — > M+ with ^i(l) = such that Gx[e,g[e)) = for all < e < 1. In view of (|ll]), g[e) 
is a monotonically (strictly) decreasing function and lime_»o5'(e) = 1. 

Now, = is the unique equilibrium solution of (^) for any e > 1. For e < \ such that 
1/(A; + 1)^ < £ < 1/A;^ holds for some integer A; > 1, one can apply the implicit function theorem 
to equation 

27r 

Gj{e, := T[e, wq) - — = , 



with j = 1, . . . ,k and conclude, exactly as in the case for Gi, the existence of a unique mono- 
tone decreasing function gj : [0, 1/j^] — > M+ with gj{l/j'^) = and \imgj{e) = 1, such that 
Gj{6, gj{e)) = 0. The branch gj{e) defines a nontrivial equilibrium solution vj of @ given by 
the p-component of the orbit 7gj(£) with {gj{e),0) at x = and p{x) = at each semi-period 
X = rmr/j, m = 1, . . . ,j. Hence, if e is such that l/(/c + 1)^ < e < there exist k equilibrium 

(not identically 0) solutions of (|^), v^, . . . , f^, with having j — I zeros in the interval (0,7r). 
Notice that gj{e) = Vj{0) > 0, j = 1, . . . , k, and Vi{x) > for all < x < vr. These will be useful 
for the stability analysis. There is an additional set of k equilibrium solution . . . with 

vT(x) =p{TT/j + x).] 

When e is small enough the orbit spends most of time at right semi-plane, p{x) ~ x mod 27r/j 
and := v^{y) dy has wrinkles of parabolic shape separated by "cusps". From (P) the 

steadily propagating solutions of (|l|) has the profile of 'dj propagating with the velocity V: 

^,{x)=^,{x)-\'^. (12) 

The plane flame front solution ipQ = Const is an asymptotically stable solution if £ > 1 for 
both equations (|1]) and (0). Its local stability can be read directly from the spectrum <7{A) of the 
operator A, given by (|I]) and (|^) linearized about ipo, which consist of simple eigenvalues 

A„ = 1 — en^ and r]n = n — en^, n = 1,2 . . . , (13) 



respectively, with corresponding eigenfunction i^l/vrcosnx. So, A„ < and ?7„ < if e > 1. 
At e = 1, the trivial solution <fQ bifurcates [Q into a steadily propagating front in which, when 
extended periodically to the real line as an even function, there is one "parabolic" tip centered at 
one wall and a single "cusp" at another. The next result of [pM|| , Theorems 5.1 and 5.14, states 
that such configuration is the only one globally asymptotically stable solution of ([I|) for all £ < 1. 
Parabolic front with centered tip or wrinkled fiame fronts may be described as quasi-equilibrium 



states discussed in refs. |[MS| , |BKS| , |SW|| . By (|]), it is enough to examine the stability of the 



nontrivial equilibrium solutions of (0). 

Ifuo G S^/^ and e > 1, then the initial value problem (Qj with u{0, ■) = uq satisfies lim \\u{t, ■)\\ii2 = 
and the trivial solution is globally stable. For e < 1 so that l/{k + lY < e < l/k"^ holds, there 



exist p > such that, if \\uo — '?^||i/2 < P; ^hen lim \\u{t, ■) — vW-^,^ = for v = vf and, for any 
sequence {Mn}„>i with lim ||u„ — f ||i/2 = 0, we have sup \\u(t, •) — f ||i/2 > 5 > for all n and v 

— n—*oo ' J- " ' 



t>0 



equal to or any equilibrium v^, j = 2,...,k. Moreover, there is an open dense set U C B^^'^ 
containing the origin and such that, if uq E U and e < 1, then lim \\u{t, ■) — ^i'||^/2 ~ 

[Let us first examine the linear stability. If u{t, x) = v{x) + C{ty then equation (|) can be written 

as 

where 

L^ = L[v]^ = ee-ve + {l-v')^ 

is the linearization of the right hand side of @ about v. Acting on the space of functions ^ satisfying 
^(0) = ^(vr) = 0, L is symmetric, (ji^L^)^ = (Lrj,^)^, with respect to the inner product 

iv,Op = / vix) ^(x) p{x) dx 
Jo 

with weight p{x) = exp [—e^^ v{y) dyY As a consequence, we can apply the comparison 
theorem to establish the following criterium (see ||G1VI|] for details). 



Stability Criterium. If is the solution of 

L[v]ip = 0, 

on < X < IT, satisfying (p{0) = and ^p'{0) = 1, the largest eigenvalue X of L [v] is negative if 
ip > on (0, tt) and positive if there exist x' such that < x' < it and ^plx') = 0. 

For an equilibrium solution t> of (|^), let 

X = c{2v-v") , 

where c > is chosen so that x'(0) = 1- It follows from equilibrium equation ev" = — {1 — v') v 
that x(0) = and x > whenever f > (recall v{0) = and (1 — f ') > for all closed orbits). 
Moreover, if f > 0, an explicit calculation gives (see [|GM|| ) 

L[v]x = ~2 {vfv<0, 
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and this implies, by applying the comparison theorem once more, 



(p{x) > x{x) 

for all < a; < TT. It thus follows from the stability criterium that is a stable equilibrium solution 
for all £ < 1 (recall Vi{x) > for all < a; < vr, implying xi^) > in the same domain). To show 
V2, . ■ . ,Vk are unstable, we observe 

L[v,]v'^ = {ev';+il-v'^)v,y = 0, (14) 

by the equilibrium equation. It thus follows that vj with j > 2 is unstable in view of the equilib- 
rium criterium and the fact that vj has at least one zero in (0, tt) and the same holds for by the 
Wronskian positivity W{ip,Vj;x) > 0. It is important to note that Vj fails to be an eigenfunction 
of L [vj] with eigenvalue because it does not satisfies the boundary conditions required. Approx- 
imate eigenfunction with exponentially small (in e) eigeinvalue can, however, be constructed using 
boundary layer techniques (see e. g. |pW|| ). Equation (|I^) holds also for any stationary solution v 
of Michelson-Sivashinsky equation. 

We turn to the global stability. A Liapunov function for equation (^ can be constructed via 
the generalized Euler-Lagrange method due to Zelenyak, Lavrentiev and Vishnevskii | |ZLV|| . Let 



U{u) = I <^{u,u,j)dx 





where $(p, w) = —y^ji^e) + (1 — w) In (1 — ty) + is an appropriated "Lagrangian" . The total 
derivative of U with respect to t is obtained by the calculus of variation 

/•TT 

= - / p{ux) i^Uxx -uu^ + u)ut dx 



where p{w) = e^^/ (l — w) is a positive weight. Note U is negative in view of (|^) and f/ is a Liapunov 
functional. Due to the fact that the trajectories {u{t, ■)}t>o ^ compact set, LaSalle's invariance 
principle can be apphed to show that all solutions of (|^) in B^^'^ converge to an equilibrium solution 
as t oo. From equations (^, (^ and the two results stated before one concludes that ipo = 



is globally asymptotically stable solution of (P for e > 1. In addition, if e < 1, there is an open 
dense set U C B^^"^ of initial condition such that all solutions 0(t,a;), with initial condition on U, is 
asymptotically of the form + Vt where t^i = Vi{y) dy, the velocity V =— /2 and vi is the 
stable equilibrium solution of (|).] 



MULTI-COALESCENT POLE SOLUTIONS 



The same scenario seems to hold for equation (0), at least if one restricts to the space of 

coalescent pole solutions. According to [P?FII] , [VM1|| , if e is such that < £ < , n = 

^Tli I 1 2?^ 1 

ex 



px 

1, 2, . . ., there exist n steadily propagating fronts, (pj{t, x) = ^j{x) + Vt, with ^j{x) = I Vj{y) dy 

Jo 

and 

vj{x) = £^cot [^-Y^j 
k=i ^ ^ 

for j = 1, . . . , n, each corresponding to the coalescent j-pole solution of (0), but only the one with 
largest number of poles, is asymptotically stable. When e crosses l/{2n + 1) from above, 
becomes unstable and a new solution (pn+i bifurcates from the former solution. 

(fj was called coalescent j-pole solution because of the dynamic of poles {zk{t)}'^fj^i, induced by 
(H), tends to align them parallel to the imaginary axis. In ||TFII| , [VM1|] , was considered with 



periodic boundary conditions. With Newmann (adiabatic) boundary conditions there are two sets 
of steadily propagating fronts {^f}^_-^ distinguished by the location of their poles: (f^ and (pj has 
poles aligned at ^ez = and ?R.ez = n, respectively. 

Note that there is a relation between the number of steadily propagating fronts and of positive 
eigenvalues in the spectrum (t{A) = {Tij}^^^ of the linear operator Av = ev" + /(f), since the 
latter increases by one unit at = l/n and the former at En = 1/ (2n — 1). For equation the 
number of steady solutions and the dimension of the unstable manifold Ain of the trivial solution 
(po = const, are equal with the instability value now located at There, in contradistinction, 

(ff are the only asymptotically stable solution for all e < 1. 



Although (|T2|) holds for (pj with Vj a coalescent pole solution of 

ev:^x -vvx + I{v) = , (15) 
10 



vf does not vanish in (0, tt). One may define a coalescent j-pole solution f j'^^^ with k — 1 zeros in 

(0, tt), A; = 1, 2, . . ., by setting v\^^^(x) = —vfikx) mod vr. Note v^^"*^ solves ([T5|) provided vf solves 

the same equation with e replaced by he. As a consequence, a sets of steadily propagating fronts 

< o?;- ^ exists if e is such that < ke < holds for some n = 1,2,... (see Figure a). 

}j=i 2n + l - 2n-l ' ' v & tj; 

In total, there are 2 [2ri/ (2m — 1)] = o(n^) coalescent steady solutions for e > - — — j-, where 

[z] means the integer part of real number z. The solutions V^j^^^ with k > 2, are not stable and 

may represent the cellular profile observed experimentally provided an associate quasi-equilibrium 

solution described in ||MS| , |BKS|, ^W|] can be defined. 

1 1 

In the following, for each e such that < £ < — holds, a new family JF„ of steadily 

n + 1 n 

propagating flame front solutions, denominated multi-coalescent n-pole solutions, will be intro- 
duced. Our preliminary investigation indicates that there are at least an exponential number c"~^ 
of solutions in JF„ and all, but 2 of them, seems to be unstable by numerical computation. As a 
consequence, the invariant set )C^^^ = [J^-Mni'^), defined as the union of the unstable manifold of 
all equilibrium solutions, for equation @ differs enormously from the invariant set )Cf^^ = A^n(v'o) 
for equation Here, n indicates the number of bifurcations with respect to the trivial solution ipQ. 
In particular, /C*^*^ may have dimension exponentially more numerous than the dimension of /C^'^ 
(for comparison, see Figures 0and|^). We believe that this crucial distinction is responsible for the 
disagreement between the numerical study by Gutman-Sivashinky ||Glj|| and the exact calculation 
by Vaymblat-Matalon ||VM1|| . 

The bi-coalescent n-pole solutions v^no.n^ are indexed by {nQ,n.,^) with hq + n.,^ = n indicating 
the number of pairs of complex conjugate polesQ] in each line z = and ^ez = n. Note that 
the system of equations governing the dynamics of the poles Zj = xj + ii/j, j = 1, . . . , 2n, in the 
complex plane, given by 

E. f Zj- zi\ . yj 
cot I : 1 — z 



l%-| 



preserves the location of real part Xj since, in this case, 3fJe cot ^^~2 — ~ ^' "^^^ poles of the 
bi-coalescent solution (Pno,n^ thus satisfies 

Xj= and yj= Fj (16) 
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where 



2 j) l%l 



with 77^7 := cos {xj — xi) taking ±1 values according the poles Zj and zi are in the same or different 
line. Because of the real parte Xj of the pole gives the "cusp" position of a propagating flame 
profile, the bi-coalescent solution has its tip centered somewhere in the interior of channel. 

Thual, Frisch and Henon ||TFH|| have proven that, provided n is such that e {2n — 1) < 1, there 
exist one and only one coalescent steady solution and any solution of ([16|) with rjji = 1 for all j and 
/, tends toward this steady state as t ^ 00. The proof of these properties is based on the existence 
of a Liapunov function with negative curvature in every direction. Here, there exists a Liapunov 
function for bi-coalescent solutions 



.2n 



satisfying f/= ^F'j > whose Hessian matrix H = [Uy^y-]jj^i cannot be proven to be 

negative definite in the case of strictly bi-coalescent solution (no, n^, 7^ 0,n/2 ) since its Gersgorin 
discs may have non- vanishing intersection with the semi-plane ^ez >0. As U may have several 



local maxima and saddle points in this case (see Figure 2), any solution of equations (|l^) tends 
toward to a steady bi-coalescent state as t — 00 but uniqueness cannot be guaranteed. 

Let the poles zi, . . . , Z2n of a bi-coalescente steady solution (pno,n^ be indexed as follows: Xj = 
if j = 1, . . . , no, Xj = IT if j = no + 1, . . . , n and i/n+j = —yj for j = 1, . . . , n. The case with 
riQ = tIt, = n/2 plays special role to describe the stability of coalescent solutions. Note that, if 
Vj = yj+n/2 for j = 1, . . . , n/2, then 

Vn/2M^) = E jcot^^^^") +COt[^^^^^^H 

j = l ,,£{-1,1} ^ ^ ^ ^ 

n/2 

= 2^^ E cot [x-iriyj) 

corresponds to a coalescent n/2-pole solution with one zero in (0,7r). As e varies from 
1/ (2n — 1) to 1/ (2n + 1), a point {zj}'^^^ satisfying yj = yj+n/2 can be shown to change from 
saddle point to a global maximum of U, turning a local maximum somewhere in between. 

12 



1 12 14 



Figure 2: Bifurcation diagram for the distance Aip from cusp to tip of steadily state solutions of 
(|lD. Solid and dashed lines refer, respectively, to stable and unstable solutions. 




Figure 3: Bifurcation diagram for the distance Aip from cusp to tip of steadily state solutions 
of (0). Solid and dashed lines refer, respectively, to stable and unstable coalescente solutions. 
Doted-dashed lines refer to unstable multi-coalescente solutions. 

Two conclusions can immediately be extracted from these observations. Because ([TBI) tends to 
align the poles along either the line 3?e z = or 3?e 2; = tt, the coalescent n-pole solution (p^ is more 
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stable than y9^^2 since (y^^^g is equivalent to a bi-coalescent n-pole solution ^n/2,n/2 which is unstable 
under small perturbation which involves the real part Xj of the poles. Moreover, we may construct 
from a coalescent n-pole solution ip^ with k — 1 zeros in (0, tt) a family of bi-coalescent solutions 
'{'no,n^/k "^ith tlq poles aligned in Sftez = and n — uq poles aligned in ^ez = ir/k which agree 
with the coalescent solution at e > 1/ {2kn + k) if hq = riT^/k = n/2 and yj = yj+n/2- Proceeding 
in a similar fashion, one can introduce muti-coalescent n-pole solutions fno,n2„/k,---n2(k-i)-K/k' with 
n2jTT/k ^ 1 poles at^ez = 2j7i/k and f^2j-K/k = ^ and which vanishes at A; — 1 points in (0, vr). 

Hence, a trajectory {zjit)}^^^ of (|16D in the phase space C", with {zj{0)} close to the poles of a 
multi-coalescent steady solution may go along many intermediate steady states before it reaches 
the final equilibrium. 



CONCLUSIONS 



Whether the cellular structure (flame profile with many wrinkles), observed numerically by 
Gutman-Sivashinky [pS|] and experimentally by [0, could be produced by the instability of the 
coalescent pole solutions for e small, has been debated in the literature (see e.g.[ [KOP| , |RAS|1 ) . 



The work of Vaymblat-Matalon ||VM1| , |VM2|| has resolved the controversies by proving that there 



always exist a unique (linearly) stable coalescent pole solution for e < 1. In [|VM1|| , the discrepancy 
between the numerical and the exact results is explained as an artefact of truncation and we shall 
not observe different profiles if more modes were included. According to Joulin [^f, once equation 
(ID is incapable to describe observed wrinkled propagating flames (as an equilibrium solution) it 
should be replaced by another model. 

A different scenario has been presented for Rakib-Sivashinsky equation. Numerical integrations 
of (p have agreed with the analytic prediction since the beginning, although truncated equation 
has been used. Besides, parabolic profile with centered tip and cellular profiles can be successfully 
explained as metastable states [[MiS| , |BKS| , |SW|| . So, the question to be addressed is why equation 
(H) is more sensitive than (|ID to be treated numerically and whether the cellular profiles can be 
described as a quasi-equilibrium solution of (H). 

Based in the present analysis of equation (|I]) and in the existence of multi-coalescent steady 
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states of (0), we argue in the following that many questions remain to be investigated before 
Michelson-Sivashinsky equation is abandoned. 

Using the analysis in [GM|, global existence in a Sobolev space Hi (for all times t > 0) and 
uniqueness can be established for equation (^). A basic question is: Does the solution (p of ^ 
with initial condition 0o in a dense subset of Hi converge, as t —>■ oo and for all e < 1, to one of 
coalescent steady state described in j\VMl^ 

To answer this question a geometric analysis, as given for equation (||), must be carried out 
for equation @. A family JF„ of steady multi-coalescent poles solutions have been described. It 
may be difficult to determine whether JF„ exhausts all steady solutions of (|^) but it is already a 
remarkable difference between both equations. If all states of except the states ip^ described 
by [[VM1|] , are shown to be unstable, then the unstable manifold JC^^ containing the origin would 
have dimension at least exponentially large with n. 

The method of [[ZLV|| may be useful to construct a Liapunov function. To establish the existence 
of a dense subspace of initial conditions from which the solutions of (|]) converge to a steady solution 
requires, besides a Liapunov function, that the trajectories remain in a compact set. For this, it 
is enough that the trajectories remain bounded, which follows if the maximum principle can be 
shown to be adapted for equation (||). 

Finally, the existence of multi-coalescent pole solutions would explain the discrepancy between 
the numerical integration by Sivashinky and the linear stability analysis in ||VM1|| . For Rakib- 
Sivashinsky equation, there exist metastable solutions whose time interval they remain "stable" 
becomes exponentially long when e is small, creating the illusion that they have reached the equi- 
librium. If metastable states can be constructed from the bi-coalescent states is a question to be 
investigate. It would, in particular, describe the quasi-stable behavior of parabolic steadily propa- 
gating flame with centered tip. Moreover, the effect of truncation would become more sensitive than 
for equation (||) in view of the fact there is exponentially more numerous (meta)states available. 
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/dy 
(j)Jt,y) is the composition of minus the Hilbert transform and the 
x-y 

derivative of cf) with respect to x. 



[2] In ||VM1|| , the Darrieus-Landau instabihty term 1(0) of (0) has a 1/2 factor in front so, the 
bifurcation point is at £ = 1/7= 1/2. In ||BKS|| , equation (|lD is defined for x G (0, 1) and 
the bifurcation occurs at e = vr"^. 



[3] Here, L is defined to be minus the L operator in ||GM 



[4] Because we are looking for real valued solutions of (0), poles must occur in pairs of complex 
conjugate. 
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